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Abstract. 

We report measurements of 57 Fe Mossbauer spectra for 7?Fe2Zn 2 o (R = Lu, Yb, 
Gd) from ~ 4.5 K to room temperature. The obtained isomer shift values are very 
similar for all three compounds, their temperature dependence was analyzed within 
the Debye model and resulted in an estimate of the Debye temperatures of 450-500 
K. The values of quadrupole splitting at room temperature change with the cubic 
lattice constant a in a linear fashion. For GdFe 2 Zn 2 o, ferromagnetic order is seen as 
an appearance of a sextet in the spectra. The 5 'Fe site hyperfine field for T —» 0 was 
evaluated to be ~ 2.4 T. 
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The AT 2 Zn 20 {R = rare earth, T = transition metal) large family of dilute, rare earth 
containing intermetallic compounds was discovered almost two decades ago, [lj but 
only a decade later Eli the physical properties of the members of this family started 
to be of interest for the broader physics community. The studies of this family brought 
to light six new Yb-based heavy fermion compounds [3], nearly ferromagnetic Fermi 
liquid behavior in YFe 2 Zn 2 o and LuFe 2 Zn 2 o Em and physics of dilute, local magnetic 
moments in nearly ferromagnetic Fermi liquid. El El For example the difference in the 
ordering temperatures and the nature of the magnetic ordered state in the GdT 2 Zn 20 
(T — Fe, Ru, Os, Co, Rh, and Ir) was understood in the model of Heisenberg moments 
embedded in a nearly ferromagnetic, highly polarizable Fermi liquid, combined with 
the electronic band filling that strongly affect this polarizability. [6] Similar approach, 
with an addition of the effects of the crystal magnetic field, can be used to understand 
enhanced Curie temperatures in the AFe 2 Zn 20 (R, = rare earth) series. [0], 7, El Ej 

The 7?r 2 Zn 20 family was mostly studied using bulk measurements of physical 
properties. Only few publications 0 nm nn reported results of the measurements 
by local probe techniques. One of such techniques is Mossbauer spectroscopy that 
allows for the determination of local electronic and magnetic properties. [12J Of the 
i?T 2 Zn 20 family, 5 'Fe Mossbauer spectroscopy was recently used to study DyFe 2 Zn 20 
and YFe 2 Zn 20 . HU The data presented for YFe 2 Zn 20 at 299 K and 78 K indicated a 
paramagnetic state of this compound, for DyFe 2 Zn 20 the data at 8 K demonstrated 
presence of the hyperfine field consistent with a long range magnetic order, El SUSIE] 
however a complex combination one one unknown subspectrum and two main subspectra 
were used to fit the data. HU 

In this work we use 5 'Fe Mossbauer spectroscopy to study the properties of 
GdFe 2 Zn 2 o, the compound that has the highest Curie temperature in the RTo Zn 2 o family 
El SI E] and has no possible complications of crystal electric field effects since Gd 3+ 
has Hund’s rule ground state multiplet J with zero angular momentum. Additionally, 
we perform measurements on paramagnetic LuFe 2 Zn 2 o, to compare with the published 
data on YFe 2 Zn 20 and evaluate the steric effects on the hyperfine parameters, and on 
the heavy fermion member of the family: YbFe 2 Zn 2 o- 

2. Experimental methods 

Single crystalline samples of RFe 2 Zn 20 (R = Lu, Yb, Gd) were grown out of excess 
Zn using standard solution growth techniques. [TH] Initial ratios of starting elements 
(R:Fe:Zn) were 2:4:94. The constituent elements were placed in an alumina crucible, 
sealed in a quartz ampule under ~ 1/3 atmosphere Ar, heated in a box furnace to 900- 
1150° C and then slowly cooled to 600° C over ~ 100 h. More details on single crystal 
growth of these materials are provided in Refs. Eli- 

Mossbauer spectroscopy measurements were performed using a SEE Co. 
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conventional constant acceleration type spectrometer in transmission geometry with an 
57 Co(Rh) source, which had an initial (8-9 months before the measurements) intensity 
50 mCi, kept at room temperature. For the absorber, the i?Fe 2 Zn 20 {R = Lu, Yb, Gd) 
single crystals were powdered using mortar and pestle and mixed with a ZG grade BN 
powder to ensure homogeneity. The absorber holder comprised two nested white Dclrin 
cups. The absorber holder was locked in a thermal contact with a copper block with 
a temperature sensor and a heater, and aligned with the 7 - source and detector. The 
absorber was cooled to a desired temperature using a Janis model SHI-850-5 closed 
cycle refrigerator (with vibration damping). The driver velocity was calibrated by ct-Fe 
foil and all isomer shifts (IS) are quoted relative to the a-Fe foil at room temperature. 
A commercial software package MossWinn o were used to analyze the Mossbauer 
spectra in this work. 

3. Results and discussion 

3.1. L'uFe 2 Zn 2 o and Y&Fe 2 Zn 2 o 

The Fe atoms occupy a single 16 d site with 3m (trigonal) point symmetry in the cubic 
CeCr 2 Al 2 o type structure (Fd3m, Z = 8 ). |T] The non-spherical environment of the 
Fe atoms (Fig. [l| ensures a nonzero electric held gradient tensor (EFG) at the 5 'Fe 
site. Indeed, the 5 'Fe Mossbauer spectra of LuFe 2 Zn 20 and YbFe 2 Zn 20 taken at different 
temperatures (Fig. [2| show similar spectral shapes with a clear quadrupole splitting. 
Each of these spectra can be fitted with a single doublet, consistent with Fe occupying 
a single, unique crystallographic site. The spectra are very similar for both compounds 
and are very weakly temperature-dependent. 

The temperature dependent hyperhne parameters for LuFe 2 Zn 20 and YbFe 2 Zn 20 
are shown in Fig. [3j These parameters are very similar for both compounds. The 
linewidth is basically temperature independent. The isomer shift increases with decrease 
in temperature. The main contribution to the temperature dependence of the isomer 
shift comes from the second order Doppler shift and is usually described by the Debye 
model: 



where c is the velocity of light, M is the mass of the 5 'Fe nucleus, and IS{ 0) is is 


the temperature-independent part, i.e. the chemical shift. The obtained values of the 
Debye temperature, are 524 ± 47 K and 487 ± 17 K for LuFe 2 Zn 2 o and YbFe 2 Zn 2 o 
respectively. 

The quadrupole splitting increases slightly, by ~ 1.5 — 3% on cooling from room 
temperature to the base temperature. 

Temperature dependence of the relative spectral area (SA) in absence of phase 
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where / is the recoilless fraction, which is proportional to the spectral area for thin 
sample and E 7 is the 7 -ray energy. The values of 0 o obtained from the fits are 410 ± 32 
K and 433 ± 25 K for LuFe 2 Zn 20 and YbFe 2 Zn 20 respectively. 

3.2. GdFe 2 Zn2o 

The 57 Fe Mossbauer spectra of GdFe 2 Zn 2 o, taken at different temperatures, are presented 
in Fig. [4j Whereas at high temperatures the observed spectra are doublets and 
are very similar to those measured for LuFe 2 Zn 20 and YbFe 2 Zn 20 (Fig. | 2 j) , at the 
base temperature the spectrum splits in six lines indicating the presence of magnetic 
hyperfine field at the Fe site, in agreement with ferromagnetic ground state of GdFe 2 Zn 2 o 
proposed based on the bulk magnetic measurements. [2] The low temperature spectra 
for GdFe 2 Zn 20 depend on the electric field gradient and magnetic hyperfine field, Bhf, 
at the 57 Fe nuclei, as well as on the angle between the principal axis of the EFG tensor 
and the direction of the Bhf. The ’’mixed Q + M static Hamiltonian (powder)” model 
(the z-axis is parallel to the principal axis of the EFG ternsor, the static Hamiltonian 
includes magnetic and quadrupole interactions with an arbitrary relative orientation) 
of the MossWinn US program was used to analyze these spectra. Due to the relative 
symmetry of the environment of the Fe atom in the lattice (see Fig. [Tj) , the EFG 
assymetry parameter, 7 = ( V xx — V yy )/V zz and the azimuthal angle of the hyperfine 
field, a were fixed at zero for the fits. 

Hyperfine field at the Fe site, Bhf as a function of temperature is shown in Fig. [5] 
At the base temperature the value of Bhf is close to 2.4 T. It has been observed |[THj 
that temperature dependence of the hyperfine field follows closely (even if not exactly) 
the saturation magnetization curve. We can fit Bhf(T ) data with the function based on 
the Curie-Bloch equation [16] for saturation magnetization. 



(3) 


where a is an empirical constant and (3 is a critical exponent. The value /3 = 1/3 
|16j was used in the fit. As the result of the fit the following values were obtained: 
Bhf{ 0) = 2.36 ± 0.03 T, Tc = 84.6 ± 0.8 K, a = 1.4 ± 0.1. The bulk measurements of 
GdFe 2 Zn 20 [E, - 6 ] yield a Tc = 86 ±2 K value of the Curie temperature, that is consistent 
with the results of the fit and with our observation of a magnetic sextet in the 80 K 
spectrum and a non-magnetic doublet in the 90 K spectrum (Fig. [4]). 

The values of the V zz component of the EFG tensor and of the polar angle f3 between 
the directions of V zz and B h f (Fig. [5]) show only insignificant change between ~ 4.5 K 


and 80 K. 
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The temperature dependent hyperfine parameters of GdFe 2 Zii 2 o are shown in Fig. 
[6} The values of IS(T), QS(T ) and the general behavior of SA(T) are consistent with 
those found for LuFe 2 Zii 2 o and YbFe 2 Zii 2 o- The overall behavior of the isomer shift 
can be described by the Debye model with ©d = 465 ± 23 K, similarly, above the 
magnetic transition, the Debye model with ©d = 450 ± 14 K can be used to describe 
the temperature dependent relative spectral area. It should be noted that an increase 
in the spectral area (beyond to what is expected in the Debye model) is apparent on 
decrease of temperature below the transition. Since significant change in phonon spectra 
is not anticipated in the ferromagnetic state, and since the IS(T) behavior in GdFe 2 Zn 2 o 
is very similar to those in LuFe 2 Zn 2 o and YbFe 2 Zn 2 o, probably some unknown artifacts 
contribute to the curious SA(T ) behavior shown in Fig. [6] The quadrupole splitting 
is only weakly temperature dependent both above and below the magnetic transition. 
QS'(295K) = 0.637±0.005 nnn/s that is ~ 10% lower than the room temperature values 
for LuFe 2 Zn 2 o and YbFe 2 Zn 2 o (0.712 ± 0.002 and 0.708 ± 0.003, respectively). 

3.3. Discussion 

The values and the temperature dependence of the isomer shift for all three compounds 
studied in this work are very similar and are consistent with the results reported for 
YFe 2 Zii 2 o and DyFe 2 Zn 2 o- m The Debye temperature, evaluated from the IS(T) data 
is ~ 450 — 500 K for all of them. Similar but slightly lower values of Debye temperature, 
~ 410 — 450 K, were estimated for RFeo^^o (R = Lu, Yb, Gd) from the temperature 
dependence of the relative spectral area. This discrepancy may be explained by the 
fact the area reflects the average mean-square displacements, whereas IS is related 
to the mean-square velocity of the Mssbauer atom. Both quantities may respond in 
different ways to the lattice anharmonicities. It has to be mentioned that the ©o 
values inferred from the temperature dependence of the hyperfine parameters above 
are significantly higher than ©£> 340 K evaluated for YFe 2 Zn 2 o and LuFe 2 Zn 2 o from 

the low temperature specific heat capacity. [0j| This discrepancy might be related to 
the fact that the ©d estimates from the temperature dependent Mossbauer spectra 
reflect, to a large extent, the phonon frequencies related to the Mossbauer atom ( 5 'Fe), 
whereas specific heat capacity involves more uniform average. Additionally, other, e.g. 
electronic, contributions might play a role in the temperature evolution of the IS and 
SA in this family of materials with nearly ferromagnetic Fermi liquid matrix. 

The 5 'Fe hyperfine parameters of the Yb-based heavy fermion YbFe 2 Zii 2 o and 
nearly ferromagnetic Fermi liquid compound LuFe 2 Zii 2 o are virtually indistinguishable. 
Possibly the effects of additional electron correlations below the Kondo temperature of 
~ 33 K in YbFe 2 Zn 2 o, [3] as probed on the Fe site, are small in comparison with the 
effects of strongly correlated Fermi liquid governed by the transition metal density of 
states. [2j 

The quadrupole splitting values are governed by the local environment of the Fe 
atoms (Fig. [Tj). To follow the geometrical changes in this environment with the change of 
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rare earth ionic radius or the lattice parameters is a cumbersome task and does require a 
detailed knowledge of atomic coordinates for each compound. Experimentally though, 
it appears that the QS at room temperature decreases as a function of the a-lattice 
parameter in close to linear fashion (Fig. [T]) suggesting the trend of the Fe environment 
to become more isotropic with increase of a. This observation is consistent with a very 
simple notion the average electric field on the Fe site becoming ’’more spherical” with 
point charges of the local environment moving further apart. 

The ferromagnetic transition in GdFe 2 Zii 2 o is seen in 5 'Fe Mossbauer spectra as 
appearance of a sextet (Fig. [4]), additionally, there is an increase in the relative spectral 
area near Curie temperature (Fig. [6]). The value of the hyperfine field for T —> 0 is 
Bhf( 0) = 2.36 ± 0.03 T, that is ~ 25% higher than ~ 1.84 T measured in DyFe 2 Zn 2 o- 
m It is possible that the crystal electric field effect present in DyFe 2 Zn 2 o and absent 
in GdFe 2 Zn 2 o should be considered in the analysis of this difference. 

The temperature dependence of the 5 'Fe site hyperfine field in GdFe 2 Zii 2 o allows 
us to address the question of interpretation of the strong deviation of the temperature 
dependence of magnetic susceptibility from the Curie-Weiss law between Tc and ~ 250 
K, discussed in Ref. f5j. Absence of hyperfine field on the 5 'Fe site at and above 90 
K, argues against the picture of formation of magnetic droplets consisting of the Gd +3 
local moments and oppositely polarized electron cloud, leaving the hypothesis of the 
temperature dependent coupling between Gd 3+ local moments as more viable. 

4. Summary 

To summarize, we performed measurements of & 'Fe Mossbauer spectra for AFe 2 Zii 2 o (R 
= Lu, Yb, Gd) from ~ 4.5 K to room temperature. The isomer shift values are very 
similar for all three compounds, and give an estimate of 0^ in the range of 480-500 K. 
The values of quadrupole splitting at room temperature change with the lattice constant 
in a linear fashion. Ferromagnetic order is seen as a sextet in the spectra. The Bhf(T) 
data can be fitted with the function based on the Curie-Bloch equation for saturation 
magnetization that yields the value of B^f at T —>■ 0 of ~ 2.4 T. No features in the 
Mossbauer spectra were associated with the Yb-based heavy fermion behavior in the 
YbFe 2 Zii 2 o at low temperatures. 

This study suggests that further work on understanding of the observed trend in 
QS within this family as well as on elucidation of the relation of the measured hyperfine 
field and bulk magnetic properties of the compounds in the i?Fe 2 Zn 2 o family is desired. 
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Figure 1 . (Color online) Local environment of the Fe atoms in f?Fe 2 Zn 2 o- (After Ref. 

my 







5 ' Fe Mossbauer effect in RFe 2 Zn 2 o 


9 



velocity (mm/s) 














5 ' Fe Mossbauer effect in RFe 2 Zn 2 o 


10 



velocity (mm/s) 


Figure 2. (Color online) 57 Fe Mossbauer spectra of (a) LuFe 2 Zn 2 o and (b) YbFe 2 Zn 20 
at selected temperatures. Circles - experimental data, lines - fits. 
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Figure 3. (Color online) Temperature dependence of the hyperfine parameters 
obtained from 57 Fe Mossbauer spectra of LuFe 2 Zn 2 o (filled circles) and YbFe 2 Zn 20 
(half filled diamonds). LW - linewidth, IS - isomer shift, QS - quadrupole splitting, SA 
- normalized spectral area. Dashed lines - Debye fits as described in the text. 
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Figure 4. (Color online) 57 Fe Mossbauer spectra of GdFe 2 Zn 2 o at selected 
temperatures. Circles - experimental data, lines - fits. 
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Figure 5. (Color online) Top panel: magnetic hyperfine field, Bhf{T), for GdFe 2 Zn 2 o- 
Symbols - data, dashed line - fit to Curie - Bloch equation - see the text. Middle 
panel: V zz componen of the telectric field gradient tensor as a function of temperature. 
Bottom panel: temperature dependence of the angle between V zz and the direction of 
the B hf . 
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Figure 6. (Color online) Temperature dependence of the hyperfine parameters 
obtained from 57 Fe Mossbauer spectra of GdFe 2 Zn 2 o. IS - isomer shift, QS - quadrupole 
splitting, SA - normalized spectral area. Dashed lines - Debye fits as described in the 
text. For IS(T ) and QS(T ) data for LuFe 2 Zn 2 o (+) and YbFe 2 Zn 2 o (x) are plotted 
for comparison. 
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Figure 7. (Color online) Quadrupole splitting at room temperature as a function of 
the a lattice parameter. Filled circles - this work, half-filled circles - from Ref. El- 
Dashed is a guide to the eye. Values of the lattice parameter a were taken from Ref. 
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